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FOREWORD 


This  document  has  been  prepared  in  order  to  provide  an  overall 
summary  of  the  work  accomplished  between  November  1974  and  June  1978 
under  Grant  No.  AFOSR  75-2784.  While  the  report  is  intended  to  be 
self-contained,  an  attempt  has  been  made  to  maintain  its  brevity  not 
only  through  appropriate  reference  to  earlier  Interim  Reports  but  also 
through  reference  to  material  arising  from  this  work  already  published 
in  the  literature.  Prof.  Sir  William  Hawthorne  of  the  University  of 
Cambridge,  England,  has  been  co-investigator  during  the  Grant  period; 
Prof.  J.E.  McCune  has  been  the  principal  investigator. 
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ABSTRACT 

A brief  account  of  the  work  accomplished  over  the  three-and-one- 
half  year  extent  of  the  subject  Grant  is  provided.  Although  the  report 
is  arranged  in  part  chronologically,  emphasis  is  given  to  that  part  of  the 
work  which  appears  not  only  to  be  new  but  also  to  have  significant  im- 
plications for  future  gas  turbine  technology.  Wherever  possible,  examples 
of  experimental  and/or  numerical  verification  of  theoretical  results  are 
noted;  shortcomings  are  also  pointed  out.  Physical  explanations  of  pheno- 
mena of  particular  importance,  of  which  some  seem  to  have  gone  previously 
unnoticed  in  turbomachinery  practice  are  offered.  Finally,  suggestions 
are  made  for  future  work  aimed  at  practical  applications  and  extensions 
of  this  research. 


I.  INTRODUCTION 

Early  analytical  studies  of  three-dimensional  flows  in  turbomachinery 
were  limited  primarily  to  what  we  shall  call  here  "strictly  linear"  theory, 
with  the  blades  of  each  cascade  being  represented  by  small  perturbations 
of  the  inlet  conditions  - usually  modeled  as  purely  axial  annular  flows. ^ 

These  initial  studies  represent  the  analogue  (for  axial  compressor  blade 
rows),  of  classical  thin  airfoil  theory.  They  were,  in  that  sense,  necessar- 
ily limited  to  very  light  wheel  loadings:  small  swirls,  small  pressure 
ratios  across  any  stage,  etc.  Such  restrictions  precluded  direct  application 
of  these  preliminary  three-dimensional  studies  to  practical  turbomachinery. 
Nevertheless,  their  results  could  not  be  entirely  overlooked,  since  they 
indicated  the  potential  existence  of  physical  effects  of  qualitative  importance 
which  had  been  up  to  that  time  omitted  in  turbomachine  analysis  and  design. 
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These  three-dimensional  phenomena  appeared  to  have  their  maximum 

significance,  according  to  the  then  available  analyses,  in  the  transonic 
1 2 

regime  ’ and  for  blade  rows  operating  with  significant  spanwise  variation 
3 4 1 

in  loading.  * ’ It  remained,  however,  to  establish  the  relative  practical 
importance  of  such  conclusions  for  turbomachinery  operating  under  realistic 
conditions,  i.e.  under  the  high  loadings  typical  of  useful  machines. 

By  the  time  the  work  described  in  this  report  was  ready  to  begin, 
the  notion  had  evolved that  it  might  be  possible  to  use  as  a reference, 
in  place  of  the  undisturbed  inlet  flow,  the  azimuthally-average  flow, 
treating  the  latter  exactly  (nonlinear ly)  with  the  help. of  the  axisymmetric 
through-flow  equations,  such  as  those  described  by  Marble.  The  0-varia- 
tion  of  the  flow  field,  it  was  argued,  could  then  be  treated  as  a perturba- 
tion about  this  "exact"  mean  flow.  If  such  a scheme  could  be  worked  out 
along  these  lines,  the  overall  large  turning  of  the  flow  would  be  described 
two-dimensionally  by  the  axisymmetric  theory  and  thus  could  be  chosen  to  be 
unrestricted  in  magnitude  through  any  given  blade  row  - pressure  ratios 
of  2.0  or  more  being  typical.  On  perturbation  of  the  axisymmetric  model, 
the  three-dimensional  structure  of  .the  flow  through  practical  (highly-loaded) 
cascades  could  then  be  described  with  sufficient  accuracy  for  a number  of 
direct  applications  within  gas  turbine  technology. 

The  basic  justification  of  this  idea  is  that  a blade  row  with  a large 

« 

number  of  blades,  B,  can  have  a large  effect  on  the  flow,  as  a result  of 
the  concerted  effect  of  all  the  blades,  whereas  each  individual  blade  is 
required  to  contribute,  by  itself,  comparatively  little. 

This  approach  to  the  analytic  study  of  axial  turbomachine  fluid  mechan- 
ics has  formed  the  basis  of  virtually  all  of  the  work  performed  under  the 
subject  Grant.  The  somewhat  intricate  methods  of  describing  effectively 
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the  perturbations  about  such  a highly  nonlinear  "reference  flow"  as  that 

» * 

mentioned  above  have  been  devised  during  this  period.  Further,  they  have 
been  gradually  improved  by  applying  the  proposed  procedure  over  a variety 
of  practical  examples,  a number  of  which  are  described  in  the  following 
sections. 

By  and  large,  these  initial  efforts  seem  to  have  been  successful  in 
providing  increased  insight  into 'the  nature  of  complex  flows  in  turbo- 
machinery. Functionally  useful  analytical  methods  have  been  developed 
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within  this  framework.  More  importantly,  their  various  applications  have 
served  to  reveal  the  importance  of  certain  physical  effects  not  formerly 
accessible  to  ready  evaluation  and/or  interpretation,  and,  in  some  cases, 
simply  unrecognized. 

As  stated  above,  the  new  theory  has  lent  itself  in  several  instances 

7 8 9 

to  an  enhanced  understanding  of  new  experimental  results.  * * In  other 

cases  it  has  stimulated  the  search,  experimentally  or  numerically,  for 

verification  of  the  predictions  it  has  provided.  The  significance  of  the 

phenomena  gradually  being  uncovered  in  this  way  is  of  potential  importance 

to  the  design  and  understanding  of  the  performance  of  practical  turbomachin- 

10,11,12 

ery  components. 

This  report  is  organized  as  follows:  In  the  following  section,  a 
brief  chronology  of  the  work  performed  is  provided.  Some  of  the  highlights 
of  the  results,  as  they  have  developed,  are  indicated  there.  In  Section  III 
certain  aspects  of  these  highlights,  and  their  possible  implications,  are 
outlined  in  more  detail. 

Finally,  in  Section  IV,  a selected  description  of  the  many  lines  open 
to  future  research  as  a result  of  the  work  accomplished  so  far  is  provided. 

An  example  is  offered  which  illustrates  some  features  of  a class  of  important 
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problems  which  can  be  treated  along  the  lines  developed  so  far.  In 
particular,  areas  in  which  the  results  can  be  extended  and  applied  in 
turbomachine  practice  are  outlined.  Some  of  this  work  has  already  begun. 
(In  fact,  at  the  time  of  this  final  writing  certain  aspects  appear  to  be 
encouragingly  near  completion.) 


II.  REVIEW  OF  PROGRESS 

The  concept  of  describing  three-dimensional  flows  in  highly-loaded 
turbomachines  by  the  method  suggested  in  the  Introduction  was  tried  first 
in  the  work  leading  to  Ref.  13.  In  that  Initial  study  an  isolated,*  moving, 
rectilinear  cascade  of  lifting  lines  with  incompressible  flow  moving 
through  it  was  used  as  a model  of  a rotor.  The  overall  turning  of  the  flow 
caused  by  the  cascade  was  taken  to  be  large,  with  the  restriction  that 
the  spanwise  variation  of  the  loading  be  small.  The  "axisymmetric"  flow 
was  replaced  in  that  model  by  the  pitchwise  average  of  the  rectilinear 
flow.  In  that  sense,  the  reference  flow  became  that  associated  with 
corresponding  highly- loaded  rectilinear  "actuator-disc"  models,  with 
weak  spanwise  variation  of  the  specified  mean  turning. 

Useful  methods  for  describing  the  three-dimensional  flow  as  a per- 
turbation of  the  above  nonlinear  (but  two-dimensional)  "reference  flow"  were 
developed  as  described  in  Ref.  13.  In  rapid  succession  the  results  were 
then  extended  (again,  for  small  variations  of  loading,  but  high  work  con- 

14 

figurations)  to  an  isentropic,  compressible  case  in  rectilinear  geometry, 
and  to  the  incompressible  annular  case.1^  The  latter  study  provided  the 
first  opportunity  for  including  important  centrifugal  effects  due  to  swirl16 
(cf.  Refs.  8,  9,11).  Somewhat  later,  the  compressible  annular  case  was 


li 


*The  assumption  of  an  isolated  row  allowed  use  of  the  approximation  that 
the  flow  was  steady  in  rotor  coordinates. 


1 


J1 

t « 

J; 

I : 


-5- 


,10 


treated  and  applied  to  a highly-loaded  transonic  cascade  of  lifting 
lines  representing  a ducted  fan  operating  near  constant  work.  These  re- 
sults confirmed  that  the  three-dimensional  effects  identified  in  the  early 
works,  i.e.  for  light  loading,*  ^ were  indeed  present  in  more  practical 
examples.  In  fact,  deviations  of  predicted  flow  angles  at  the  blade  row, 
relative  to  those  expected  in  the  absence  of  significant  three-dimensional 
effects,  were  magnified  by  increased  loading  (Fig.  1),  roughly  in  accordance 
with  expectations  arising  from  strictly  linear  theory. ^****®’*^ 

After  the  extension  of  the  original  work  to  the  annular  case  (incompressible. 
Ref.  15)  it  became  possible  to  compare  our  predictions  with  data  obtained  from 
experimental  study  of  flow  behind  a 23-bladed  transonic  rotor. ^ The  tests 
referred  to  there  were  carried  out  in  the  MIT  Blowdown  Compressor  Facility.*^ 

(The  comparison  between  a theory  of  incompressible  flow  behind  an  annular 
cascade  and  data  from  a transonic  blade  row  were  made  possible  by  adjusting 
the  annular  area  in  the  former  case  to  obtain  the  proper  mean  axial-velocity 
ratio  corresponding  to  the  mean  density  and  area  changes  that  occurred  in  the 
actual  experiment.)  Since  we  were  primarily  interested  in  flow  effects  induced 
by  vorticity  present  in  the  blade  wakes  downstream  of  the  rotor,  compressibility 
effects  other  than  those  grossly  affecting  the  mean  flow  were  not  expected 
to  be  of  primary  importance.  (As  will  be  seen  subsequently,  this  assumption 

e 

is  not  entirely  justified,  since  the  interaction  between  pressure  and  vorticity 
disturbances  is  important  in  swirling  flow.) 

Results^  of  this  comparison  were  sufficiently  encouraging  to  warrant 
continued  research  in  the  established  framework.  At  the  same  time,  not 
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These  results  led  to  Increased  efforts  toward  obtaining  a more  complete 
theoretical  description  of  actual  flows  through  practical  axial  turbo- 
machinery. 

Meanwhile,  the  work  accomplished  up  to  that  point  was  reported  in 

18 

the  first  Interim  Report  under  the  Grant.  At  about  the  same  time,  our 

results  were  summarized  in  an  invited  paper  at  a Project  Squid  (-NASA-AFOSR) 

19 

workshop,  and  later  published. 

It  should  be  emphasized  once  more  that  the  three-dimensional  analyses 
reported  in  Refs.  18  and  19  correspond  to  airfoil  and/or  propeller  theory 
in  an  ideal  gas,  with  the  important  refinement  that  the  perturbations 
introduced  by  the  individual  blades  are  pictured  as  being  superimposed  on 
a complex,  swirling  (axisymmetric) "ref erence  flow"  rather  than  on  a nearly 
uniform  one.  This  meant,  for  example,  that  in  the  initial  works  the  vorticity 
in  the  "wakes"  behind  the  blades  was  limited  to  that  associated  with  spanwise 
variation  of  circulation  along  the  blades.  Such  loading  variations  lead 
to  trailing  vorticity,  or  shed  circulation,  in  the  wakes,  much  as  in  classical 
airfoil  theory.  This  type  of  vorticity  is  aligned  with  the  (mean)  down- 
stream flow  and  is  often  called  "Beltrami  vorticity".  When  it  is  present 
to  a sufficient  extent  it  produces  a type  of  structure  in  the  velocity 
field*  quite  distinct  from  that  associated  with  the  lossy-viscous  wakes 
normally  considered  dominant  in  turbomachinery.  It  was  accordingly  of 
decisive  Importance,  in  the  comparisons  made  in  Ref.  7,  to  discover  whether 
or  not  such  velocity  field  structures  were  indeed  observed  in  the  flow 
'behind  a nonun if ormly-loaded  rotor. 
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*The  expected  structure  included,  for  example,  sawtooth-like  variations 
in  the  radial  velocity. 
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In  the  experiment  ’ these  unique  structures  were  found  to  be 
clearly  discernible  downstream  of  the  rotor  at  radial  stations  for  which 
there  was  a sufficiently  large  spanwlse  variation  of  loading  on  the  blading 
(Figs.  1-4).  At  the  same  time,  at  radial  stations  for  which  the  blade 
loading  variation  was  small  (Fig.  1),  the  wake  flow  structure  (Figs.  5-7) 
was  found  to  be  dominated  by  that  more  typical  of  viscous  wakes  behind 
airfoils  or  cascades. 

The  conclusions  drawn  from  this  comparison  were  that  the  theory  at 
that  stage  could  be  considered  (partially)  successful,  in  that  it  accounted 
for  a significant  portion  of  the  observed  wake  structure,  and  that  the 
experiment  thus  provided  strong  evidence  for  the  presence  of  downstream 
trailing  vorticity.  A three-dimensional  theory  of  the  type  being  developed 
was  therefore  necessary  in  order  to  contribute  to  estimates  of  the  flow  angles 
at  the  individual  blades,  since  such  estimates  now  had  to  include  the  in- 
duced angle  associated  with  the  "Beltrami"  vorticity. 

However,  the  comparison  made  in  Ref.  7 also  indicated  that  there  were  sub- 
stantial regions  downstream  of  the  rotor  where  the  wake  flow  structure  was 
found  to  be  dominated  by  effects  more  typical  of  viscous  wakes  behind  the 
blading.  This  is  describable  in  terms  of  the  corresponding  vorticity  (which, 
unlike  the  Beltami  Vorticity,  arises  because  of  viscous  interaction  between 
the  fluid  and  the  solid  surfaces  formed  by  the  blades)  which  is  swept  into 
these  wakes  from  the  flow  affected  by  the  blade  boundary  layers.**  Study  of 
this  effect  were,  in  large  part,  motivated  by  the  work  of  Kerrebrock*^’^  which 
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shows  that  vorticity-and/or  entropy-perturbations  must  couple  significantly 
with  pressure  (and  hence  flow-angle)  disturbances  in  highly  swirling  flow, 
•jnthe  resulting  aerodynamic  modelling  designed  to  include  this  effect,  the 
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swrirling  downstream  flow  was  pictured  as  being  threaded  with 

vortex  filaments  (oriented  in  the  spanwise  direction  immediately  aft 

of  each  blade)  lying  within  the  blade  wakes.  Because  the  vorticity 
field  can  be  related  to  the  variation  of  the  thermodynamic  properties  of 

the  fluid,  the  presence  of  such  viscous  wakes  necessarily  implies  a 

variation  of  entropy  from  streamline  to  streamline,  or,  in  the  incompressible 

limit,  a corresponding  variation  of  stagnation  pressure. 

For  this  reason,  the  available  blade-loss  coefficients  measured  for 

cascades  operating  at  typical  incidence  angles  were  used  to  establish  the 

strength  of  wake  vorticity  of  this  type.  The  combination  of  a strong  centrifugal ly 

accelerated  field  with  the  vorticity  arising  from  the  blade  wakes  turns  out 

to  lead  (perhaps  surprisingly,  at  first  glance)  to  significant  flow  angle 

11  12 

changes  at  the  blade,  (Fig.  8)  * . In  turbomachinery,  in  other  words, 

there  is  an  important  "induced  incidence"  associated  with  the  viscous  wakes 
themselves ! 

Meanwhile,  Greitzer^  had  developed  an  approxirtfte,  beautifully  intuitive, 

analysis  which  includes  many  of  the  important  three-dimensional  flow  features 

arising  from  the  interaction  of  a total  pressure  defect  (e.g.  due  to  an 

inlet  distortion)  with  a swirling  mean  flow.  His  results  were  verified  by 

« 

experimental  data  obtained  using  an  annular  rig  constructed  at  the  Whittle 

g 

Laboratory,  Cambridge  University,  England.  In  that  experiment  the  existense 


fc 


of  a systematic  downstream  static  pressure  variation  was  observed.  This  work 
motivated  even  further  development  of  our  theory,  based  on  the  analytical 
techniques  applied  in  Ref.  11.  That  work  directly  complements  and  extends 
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Greitzer's  prediction,  and  also  allows  more  detailed  comparison  with  the 
experimental  data  (Figs.  9,10,11,12).  For  the  lowest  circumferential  harmonic 
(i.e.  n=  1)  representing  the  total  pressure  defect,  there  is  substantial  agreement 
with  Greitzer's  analysis,  especially  in  the  region  not  too  far  downstream  from 
the  blade-row. 

Even  so,  for  free-vortex  swirl,  the  more  complete  theory*^  clearly  also 

shows  that  while  the  resulting  downstream  disturbances  (including  the  static 

pressure)  grow  initially  in  strength  for  a certain  distance  downstream  of  the 

blade-row,  they  then  decay  inversely  with  (some  power  of)  the  axial  distance 

downstream  (Figs.  13,14,15).  This  analysis  also  demonstrates  that  the  amplitude 

of  the  downstream  disturbances  decreases  as  the  number  of  significant  circumferential 

* 

harmonics  of  the  inlet  distortion  increases  (compare  Fig.  16  with  13).  It  is 
essential  to  note  that  these  fluid  mechanical  phenomena  are  beyond  the  predictive 
capability  of  the  two  dimensional  cascade  theory.  They  are  simply  absent  in  that 
description. 

The  work  under  the  Grant  accomplished  up  to  this  point  was  reported  in 

26 

the  second  Interim  Report  submitted  to  AFOSR.  Further,  results  on  the  inlet 

9 

distortion  problem  were  published  in  the  ASME  transactions,  while  initial 

11  26 

results  on  the  three-dimensional  effects  of  viscous  wakes  * were  published  in 
Ref.  12. 


Up  to  the  stage  described  above  we  had  confined  our  analysis  to  an 
Isolated  blade-row  (or  single  actuator-disc)  model  involving  not  only  steady 
flows  but  also  nearly  "free-vortex"  downstream  swirl.  However,  one  of  the 
important  results  of  the  comparisons  with  experimental  data  in  Ref.  7 was  to 


*ThJ.s  suggests,  of  course,  that  the  significance  of  the  effects  discussed  above 
may  diminish  in  relative  importance  as  the  number  of  blades  in  the  corresponding 
blade  row  is  increased. 


\ 


suggest  the  possibility  of  removal  of  some  of  these  restrictions.  It  was 


noticed,  for  example,  that  the  original  Beltrami-flow  calculations 
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were  successful  even  in  regions  departing  significantly  from  free-vortex 


swirl  (constant  work) . This  strongly  suggested  that  the  basic  procedure 


used  for  describing  the  three-dimensional  b lade- to-b lade  flow  (due  to  the 


presence  of  discrete  wake  structure  of  any  type)  relative  to  a mean  swirling 


flow  of  high  loading  was  in  no  way  limited  in  principle  to  rotors  with  small 


radial  variation  in  blade  circulation.  This  confirmed  that  the  basic  idea 


underlying  the  suggested  "perturbative  technique"  is  related  to  the  fact. 


mentioned  in  the  Introduction, that  a blade-row  of  a large  number  of  blades. 


B,  can  have  a large  effect  on  the  flow  even  though  each  individual  blade  is 
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only  lightly  loaded  comparatively.  ’ On  this  basis,  a more  ambitious 


analytic  study  of  strong  Beltrami  flow  was  undertaken  successfully.  Comparison 


of  this  newly-developed  analysis  with  the  available  experimental  data^'^ 


obtained  from  M.I.T.  blow-down  compressor  facility  was  then  carried  out.  That 


study  confirmed,  even  more  strongly,  the  trend  apparent  in  Ref.  7.  (Fig.  17  & 18). 


Finally,  the  opportunity  to  extend  the  analytical  techniques  we  had 


developed  up  to  that  point  under  the  Grant  to  the  treatement  of  unsteady  flow 


through  turbomachine  blade  rows  became  apparent.  A study  of  three-dimensional 


effects  in  unsteady  flows,  resulting,  for  example,  from  the  passage  of  distorted 


inlet  flow  through  a (rectilinear)  moving  cascade  was  undertaken.  The 


three-dimensional  "bending  and/or  torsional"  motions  of  the  blades  (in  rectilinear 


geometry)  have  now  been  analyzed  in  terms  of  the  three-dimensional  characteristics 


of  the  corresponding  unsteady  flow.  Initial  work  on  three-dimensional  flutter 


analysis  indicates  that  the  flutter  boundaries  of  a compressor  array  may  be 


changed  significantly  as  a result  of  unsteady  three-dimensional  effects. 
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III  EXAMPLE  OF  THE  GENERAL  ANALYTICAL  TECHNIQUES 
In  this  section,  we  present  an  illustration  of  our  general  techniques 
as  applied  to  the  study  of  unsteady  three-dimensional  flows  through  typical 
turbomachine  blade  rows  of  the  annular  and  rectilinear  type. 

The  Euler  momentum  equations  for  an  incompressible  fluid  can  be 
written  in  "Lamb's  form"  as: 


- =■  - v xn  - -v(~)  + f 

3t  P 


(3.1) 


where  V is  the  velocity,  p the  stagnation  pressure,  = V x V the 
(absolute)  vector  vorticity,  t the  time  and  +F  is  a body  force  per  unit  mass 
(local  imposed  acceleration)  reacting  on  the  fluid  as  a result  of  forces 
concentrated  at  blade  surfaces. 

In  the  absence  of  sources  or  sinks  the  continuity  condition  is,  for  this 
example,  simply 


V • V = 0 


(3.2) 


Let  us  assume  that  a velocity  potential,  <f>,  represents  the  irrota- 
tional  part  of  the  flow  field,  while  a vector  A describes  the  solenoidal  part. 
The  effect  of  bound  vorticity  on  the  fluid  will  be  introduced  by  the  use  of 
the  "body  force"  denoted  by  F.  As  pointed  out  above  F represents  the 
reaction  on  the  fluid  of  the  forces  experienced  by  the  blades.  Clearly,  then, 
away  from  the  blading,  F is  zero.  Thus,  it  is  useful  to  write,  with  V = V<J>  + A, 
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n - V X A *=  (Vo  X Vt)  + (Va  x VG)H'(a)  (3.3) 

18  22  26 

employing  a "multiple"  version  of  the  Clebsch-Hawthorne*  approach.  * * 

The  idea  of  using  more  than  one  combination  of  the  type  Vo  x Vt,  etc.,..., 

26 

was  first  found  to  be  useful  in  steady  flow  applications.  Perhaps  its 

greatest  value  lies  in  making  possible  the  identification  of  various 

quantities  with  the  physical  causes  for  the  presence  of  various  types  of 

22 

vorticity,  a point  emphasized  by  Farokhi. 

A Ojr 

In  (3.3)  H (a)  = — H (a)  , a = 0 - f (r,z,t),  where  H (a)  is  the 

D 

unit  "staircase  function"  with  jumps  of  +1  at  the  blade  surfaces 
2ir  4n 

a = 0,  ± — =■  , ± — , ...,;  B is  the  number  of  blades  in  the  blade  row  to  be 
IS  D 

described.  At  a = 0±,  we  have  H±  = ± 1/2,  or  H±  = i 7 . H'(ct)  is  then 

D 

(2ir/B)  6 (a),  where  6 (a)  is  the  "periodic  delta  function",  zero  except  at  the 
p P 

blades.  For  purposes  of  convenience  below  we  note  that  the  "sawtooth  function", 

A 

S(a),  (familiar  to  propeller  theory)  is  identically  (H(a)  - a)  and  hence 

2*jr 

S*(c0  = — 6 (a)  - 1.  Note,  then,  that  both  S(a)  and  S' (a)  have  zero  azimuthal 

a r 

average  around  the  annulus  and  that,  unlike  H(a),  S(a)  is  non-secular  in 

A 

0,  i.e.,  any  "cyclic  constant"  introduced  by  a function  like  H(a)  is  removed 
through  the  use  insterd  of  S = H(a)  - a.  Finally,  note  that  all  these 
(singular)  functions  are  to  be  understood  in  the  sense  of  generalized  functions. 
Thus,  for  example,  the  second  term  on  the  rhs  of  (3.3)  is  identically 
VH(o)  x VG. 


* The  Clebsch  method  (see  Ref.  24)  can  be  thought  of  as  a particular  choice 
of  "gauge"  for  the  vector  potential,  A.  Its  special  effectiveness  for 
flow  in  turbomachinery  was  emphasized  by  Hawthorne. 23 


Substitution  of  (3.3)  in  (3.1),  together  with  expansion  of  the  vector 
triple  products  in  V x Q yields,*  after  some  rearrangement: 


£vT-j£va+fvG-fvH(a) 


- V<7T  + If  + °!t  + ® !£>  + f (5.0  • 


(3.4) 


where  D/Dt  = 3/3t  + V • (3/3x).  We  note  first  that  DH/Dt  = 

(2tt/B)  5 (a)  (Da/Dt),  , and  this  term  vanishes  because  of  the  boundary 

p be  ,w 

condition  at  the  blades  (or  wakes)  (Da/Dt),  , = 0.  (6  (a)  serves  the  role 

bl,w  p 

of  "picking  out"  those  values  of  a = ct.  which  define  the  blade  (wake) 
surfaces.  It  is  on  such  surfaces  that  the  concentrated  (bound  and  "shed") 
vorticity  associated  with  the  blade  forces  and  their  time  - and  spanwise  - 
variation  appears.  (The  extent  to  which  this  "concentrated"  vorticity 
actually  becomes  diffuse  of  course  depends  on  the  relative  importance  of 
viscous  and  turbulent  effects.) 

Further,  on  physical  grounds  within  the  present  model,  we  expect  F(x,t) 

to  be  of  the  form  aVH(a)  = a(x,t)  6 (a)Va,  i.e.,  concentrated  at  the  blade 

P 

surfaces  (and  — except  for  viscous  stresses  — directed  normal  to  those 

surfaces  at  any  instant) . We  should  also  expect  F to  be  related  to  the 

« 

acceleration  of  an  identifiable  physical  quantity.  Thus,  we  can  reasonably 
choose  G and  F so  as  to  put 


F = -(DG/Dt)  VH(a). 


(3.4a) 


A 

* Note  that,  with  (3.3),  y = V<J>+  oVt  + HVG,  (plus  any^ gradient,  as  determined 
by  boundary  conditions'*.  Alternatively,  defining  <f>  = <J>-  aG,  we  can  write 
V = V4>  + oVt  - GVtt  + (H(ct)  - a)VG,  reflecting  the  required  non-cyclic 
nature  of  V. 
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As  we  shall  see  below  this  is  a key  to  finding  a useful  solution  of 
(3.4),  since  it  turns  out  to  define  "G"  in  a natural  way,  and  also  makes 
contact  with  classical  results  for  unsteady  thin  airfoil  and  linearized 
cascade  theories.* 

Interestingly,  the  interpretation  of  "G"  and  "a"  (in  F « aVH)  turns 
out  to  be  relatively  straightforward  in  turbomachinery:  -|^a  is  the 
incremental  local  specific  torque  introduced  at  each  blade,  while 
(-G)  is  the  mean  (azimuthally  - averaged)  specific  angular  momentum  of  the 
fluid  introduced  by  all  the  blades  up  to  any  axial  station. 

In  any  case,  with  the  above  use  of  the  blade  boundary  condition  and  the 
definition  (3.4a)  of  the  relation  between  F and  DG/Dt,  (3.4)  reduces  to 


V.V 


25/T  . W +~  +H  + o|t  + Hfa)  |£) 


3t  v3t 


3t 


(3.5) 


of  the  form:  Dr/Dt  = 1 , Va  = (—  + ■§£  + + h|£)  and 

p dt  dt  dt 


Dt  “ °* 


(3.6) 


t 


I 


i 

i 


I 


t 


* From  the  point  of  view  of  classical  potential  theory,  "G"  plays  the  role 
of  a cyclic  constant  which  must  be  specified  (because  of  the  annular 
geometry)  to  determine  the  flow  uniquely.  In  fact,  just  behind  a blade 
row,  (2tt/B)G  is  the  instantaneous,  individual  blade  circulation  at  the  : 

corresponding  spanwise  location.  As  will  be  pointed  out  later,  the 

requirements  needed  for  defining  a self-consistent  mean  flow  imply  an  i 

additional  constraint  on  G (beyond  the  definition  (3. 4. a))  which  is 
similar  to  a Kutta  condition. 


Actually,  one  can  show  that  this  is  the  general  solution  of  (3.4),  once 
it  has  been  reduced  to  (3.5).  Consistent  with  earlier  remarks,  we  can 
also  write  the  conserved  quantity  as 


a - £)c 


1 p at 


°§! ) 


(3.6a) 


The  conserved  quantity,  a , in  (3.6)  or  (3.6a)  is  the  generalization, 
for  unsteady  incompressible  flow  including  vorticity,  of  the  more  familiar 
Bernoulli  constant  of  potential  flows.  Note  that  o reflects  the  presence 
of  bound  vorticity  (non-zero  blade  forces)  and  its  time  variation  (varying 
blade  circulation)  through  the  (3G/3t)  term.  More  specifically,  it  reflects 
the  specific  power  put  into  or  extracted  from  the  fluid  by  the  blade  forces 

/v 

F at  the  moving  surfaces,  a(t).  (It  can  be  shown  that  F • V = -a^2-  = 

~ ~ a t 

3ot 

-a  6 (a)  XT’  for  this  case.) 

p a t 

The  requirement 


Dr 

Dt 


1 


(3.7) 


Identifies  "t"  as  the  fluid-particle  drift  time  (see,  for  example,  the 
discussion  given  by  Hawthorne,  Ref.  23).  The  general  solution  of  (3.6)  can 
be  found  with  the  help  of  T.  . 

To  illustrate,  consider  the  situation  for  which  the  unsteadiness  of  the 
flow  is  introduced  by  the  blade  row  (or  rows)  itself,  so  that  for  the  present 
example,  the  flow  is  steady  "far  upstream"  (e.g.,  in  the  inlet).  Then, 
from  (3.6)  and  the  definition  of  a 


o(x,t)  = p (x  (x,t))p  + const. 

^ O »oo  <w 


where  Xj^in  the  upstream  position  from  which  the  fluid  particle  at 
(x,t)  emanated  earlier.  Considering  the  meaning  of  T,  as  implied  by 
(3.7),  we  can  write 


• T (x , t) 

dx'  V (x\x’) 


where 


?’<*  , T’)  = X_00  + 


dx"  V(x",X") 


Thus,  while  Vx  = 0x^/3x)  • Vp""*^  (p-1)  , and  p— 


3x 


*r-o>  _ — o 

• ,po 


3*-  ® 

_ = I - VxV(x,t);  ~~  = -v|| 


so  that,  from  (3.8)  with  (3.7),  we  have  p|£  = - |l  v • Vp  "“and 

ot  dt  >v  * o 

3x_£d  w 

pV  • V°  = V . (3—  ) • Vp pv  * Vo  = (1  - V *Vx)  V • vp— 

%r  y *“  *•  O 


whence  Dc/Dt  = 0 for  the  circumstances  stated.  In  that 


7o  - r Vpo"(5— <;•*)) 


case 


_ 1 

I;  1 
1 1 


: 


1 

I 


h 
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.and  that  part  of  the  vortlclty  in  (3.3)  which  depends  on  Vo  is  absent 
whenever  the  inlet  flow  is  uniform. 

23 

These  basic  concepts  are  familiar  from  the  theory  of  secondary  flow  ; 
for  this  reason,  we  refer  to  the  (Vt  x Vo)  - term  in  (3.3)  as  "secondary 
vorticity",  and  denote  it  by 

Using  the  requirement  (3.2)  and  the  forms  of  V appropriate  to  (3.3) 

(see  footnote  following  Eq.  (3.4))  we  have 

O n A A A 

V 4>  = - Vo  • VT  - oV  T - VH  • VG  - HV  G 

or  (3.12) 


2 ? 9 

V <j>  = - Vo  • Vt  - oV  x - VS  • VG  - SV  G + Va»VG  + GV  a 

A A A 

where  $ = $ - aG  and  S(a)  = H(a)  - 1 . Since  a itself  depends  on  ({)  (or  <{>), 

Eqs.  (3.12)  are  non-linear.  They  are,  however,  in  a form  appropriate 

to  relatively  straightforward  perturbative  methods  with  regard  to  a 

sufficiently  well-chosen  "reference  flow." 

In  the  following,  for  illustrational  purposes,  we  limit  ourselves  to  the 

00 

case  of  "undistorted"  inlet  flow,  i.e.,  (Vp  ) = 0.  In  that  situation 

o 'x 

00 

the  only  vorticity  present  is  that  due  to  the  blade  forces  and  their  time  - 
and  spanwise  - variations.  These  latter  effects  lead  to  "shed"  and/or 
"trailing"  vorticity  in  the  wakes  of  the  kind  discussed,  for  example,  in 


. 


Refs.  7,  21,  20,  13  and  in  the  Introduction. 


III.  A.  UNIFORM  INLET,  INCOMPRESSIBLE  FLOW 


In  this  particular  case  Do/Dt  = 0 **  Vo;  thus,  throughout  the  flow. 


Alternatively,  when  O = const# , the  term  ovt  in  V and  the  term  o*; 
the  Lamb  definition  of  O can  simply  be  included  in  the  corresponding 
contributions  of  the  potential,  <}>.  Thus,  with  a = O = ct,  let 


whence 


Here,  o is  simply  a "renormalized  Bernoulli  constant",  reflecting  the  absence 


of  inlet  distortion  effects 
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The  formulations  discussed  above  are  exact  (for  the  stated  physical  model) 
but  are  not  particularly  amenable  to  further  useful  physical  interpretation. 
Before  turning  to  perturbation  techniques,  however,  we  can  extract  at 
least  one  additional  item  of  significant  information.  If  we  choose  the 
reference  flow  (or  "primary  flow"  — see  Ref.  23)  to  be  the  axisymmetric 
equivalent  (azimuthal  average)  of  the  actual  flow,  we  must  recognize  that 
such  a "reference  state"  cannot  be  time- independent  in  general.  For 


turbomachinery  of  practical  loadings  any  significant  time  dependence  of  the 
actual  flow  is  reflected  strongly  in  its  (axisymmetric)  azimuthal  average.*  The 
time  average  of  any  flow  quantity  is  not,  in  general,  equal  to  its  6-average, 
and  vice-versa. 

III.  A.  i.  THE  MEAN  (0 -AVERAGED)  FLOW  AS  REFERENCE 
In  the  uniform-inlet,  incompressible  case  the  mean  flow  is  defined  by 


H-!x?  = v(f)+!o  * 


(3. A. 2) 


where  ( ) =(l/2ir)  d0  ( ),  and 


V • V = 0 . 


The  mean-flow  absolute  vorticity  is 


n = Va  X VG 
-s 


(3.A.3) 


(3.A.4) 


where,  with  a = 0 - f(r,z,t)  , we  must  require  G = - rvfl.  (Since  we  expect 


* The  important  exception  to  this  statement,  of  course,  is  for  situations  in 
in  which  the  time  dependence  can  be  transformed  away  with  a simple  Galilean 
transformation.  For  the  "isolated  rotor"  problem,  for  example,  quantities 
depend  only  on  0'  = 0 - cot. 
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no  secondary  vorticity  in  this  example,  (3. A. 4)  contains,  as  a key  part  of  the 
definition  of  the  mean  flow,  simply  the  0-  average  of  the  last  term  in 
(3.3).)  The  meaning  of  Fo  will  be  discussed  below. 

Inserting  (3. A. 4)  in  (3. A. 2)  and  expanding  as  before,  we  find 

IvG-fS*,.-V(^  + tf-G|t>+7  (3.A.5) 

where  now  (D/Dt)  = 0/3t  + V • V)  and  where  we  have  used,  consistently 
with  (3. A. 4), 

v = V$  - GVa  = Vijr  + aVG.  (3. A. 6) 


Correspondingly,  in  view  of  (3. A. 3),  (cf.  (3.12)): 


v2*  - -Va  • VG  - aV2G 
or 

V2<J>  = +VG  »Va  + GV2a 


(3.A.7) 


The  0-component  of  (3. A. 5)  reveals 


<rF0>o 


dg  3(rV  _ 

5T  s aT  + Y * v (rve) 


(3.A.8) 


Thus,  recalling  the  form  of  the  mean  angular  momentum  conservation  equation, 

we  see  that  (rF0)  is  identically  the  mean  specific  torque  applied  to  the 
o O 

fluid  as  the  result  of  the  concerted  action  of  all  the  blades  in  a given 


1 

I 

I 


I 


r 


i 

j 


j 

; 


/ 


I 
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blade  row.  This  does  not  yet  imply,  however,  that  Fq  is  the  0-average  of 
F in  (3.1)  or  (3.4).  As  (3. 4. a)  shows,  in  fact,  (with  a = 0 - f(r,z,t)) 


F - - (|f  + (V  + v)  • VG)bl>(,  Vc, 


(3.A.9) 


or 


rF0  *=  - -<v>  • VG;  (G  = -rvQ)  (3. A. 9a) 

- - a, 

Here,  v is  defined  through  V = V + v,  and  w is  its  local  value 

at  each  blade  or  wake  surface.  We  then  see  that  rFQ  and  (rFQ)  differ  unless 

o o o 

<v>  • VG  = 0 

bl,w 

i.e.,  unless  the  meridional  component  of  the  variation  of  V about  its 
0-average,  T,  either  vanishes  at  each  blade  (or  wake)  surface  or  itself 
lies  in  surfaces  of  constant  mean  angular  momentum  (constant  G) . (Such  "G- 
surfaces"  are  surfaces  of  revolution  about  the  hub  of  the  annulus;  it  is 
sometimes  helpful  to  think  of  G as  a "Stokes  stream  function"  for  the  mean 
vorticity,  together  with  the  corresponding  physical  meaning  of  such  quantities. 
Thus,  G measures  the  flux  of  vorticity  through  surfaces  enclosed  by  G = const., 
and  is  then  the  circulation  for  corresponding  paths  enclosing  the  hub.) 

Consider  now  the  0-average  of  Eq.  (3.1): 


(3.A.10) 


3V  — p“ 

+ 


(3. A. 11) 
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In  the  present  example,  this  averaging  process  can  be  carried  out  without 
approximation,  even  for  the  non-linear  term,  because  of  the  concentrated 
("singular")  nature  of  fl.  The  result  is 

3V  _ P 

~ + V x ft  + <v>  x ST  = — V ( — — ) + F (3. A. 12) 

3t  ~ ~ bl,w~  P 

Comparing  this  exact*  average  with  the  definition  (3. A. 2),  we  see  that 


/ 


F - F = <v>  x = (<v>  • VG)Va  - (<v>  • Va)VG  (3. A. 13) 

~ ~°  bl,w  ~ bl,w  bl,w 

When  (3.A.10)  holds,  the  8 -components  of  Fq  and  F are  the  same,  consistent 

with  their  both  producing  the  same  mean  meridional  vorticity  [of  the  form 

1 3G  1 3G  — 

(—  - — —■)  ] and  the  same  mean  specific  angular  momentum,  rvfi  = -G. 

X dZ  y • • • 0 T dr  v 

Thus,  (3. A. 10)  must  hold  for  the  sake  of  self-consistency.  This  is 

tantamount  to  specifying  the  required  cyclic  constant  for  the  potential  part 

of  the  mean  flow,  as  mentioned  earlier. 

Then,  if  there  is  a remaining  difference  between  F and  Fq,  it  lies  in 

meridional  planes  and  (locally)  along  VG.  The  coefficient  is  (<v>  • Va) , 

bl,w 

and  this  is  related  to  the  boundary  conditions  at  each  blade  (wake)  surface. 

0 

At  such  surfaces  we  recall  that  (Da/Dt)  = 0,  i.e., 

bl,w 


1 


3a 

3t 


v 


Va  + <v>  • Va  ■ 0 

~ bl,w 


(3.A.14) 


' 4 


* "Exact"  within  the  present  physical  model. 


On  the  other  hand,  the  requirement  (3. A. 3)  allows  us  to  define,  even 


though  the  flow  may  be  unsteady,  a set  of  surfaces  (cxq , ) such  that, 
instantaneously,* 

V “ Va  x ViJ>  + cor§Q  (3. A.  15) 

~ o a 

with  a = 0 - f (r,z,t),  ij>  = ijj(r,z,t).  When  the  blade-row  solidity  is 
o o 

sufficiently  large  (i.e.,  when  a sufficiently  large  number,  B,  of  blades 
are  doing  the  required  job  of  a given  blade  row)  we  can  expect  the  "swirl" 
surfaces,  aQ,  to  coincide,  very  nearly,  with  the  surfaces,  a,  in  which  the 
vorticity  is  concentrated.  (In  the  limit  B -*■  00 , the  a and  a surfaces  must 
merge,  a notion  which  forms  the  basis  for  the  use  of  axisymmetric  through-flow 
theory.)  Let  us  then  define 

C(r,z,t)  = a - a = f (r,z,t)  - f(r,z,t)  * (3. A. 16) 

o o 

with  the  expectation  that  C -*■  0 as  B -*■  00 . 

Clearly,  from  (3. A. 15),  at  any  instant 

V • 7a  = to  (3. A.  15a) 

V O 

whence,  using  (3. A. 14),  the  blade  (wake)  boundary  condition  yields. 


* We  include  explicitly  the  divergence-free  vector,  coreg,  to  allow  the 

possibility  that  the  a-  surfaces,  and  hence  a -surfaces,  may  be  rotating  at 
angular  speed  w,  relative  to  "absolute",  i.e.,  duct-fixed,  coordinates. 
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<»«  ar  _ 

0)  + ~ + |r  + V • VC  = -<v>  • V (a  + C) 

3t  3t  - - bl,w  o 


(3. A. 16a) 


_=_<v>  • V (a)  -(rjt  + u) 

bl,w 


Now,  returning  to  (3. A. 5)  and  using  (3. A. 9),  (3. A. 13),  and  the  constraint 


(3. A. 10)  together  with  (3. A. 16),  we  can  write:  [recall,  a = - DG/Dt  (=  -DG/Dt, 


from  (3.A.10))] 


F = aVa  + ®VG 
~o  Dt 


(3.A.17) 


so  that 


or  VG  + §£  VG  - 0 


and,  hence,  since  V • Va  = w, 
~ o 


°--7'T  + i - “if* 


(3.A.18) 


Thus,  in  fact,  the  mean  "Bernoulli  pressure"  is  modified  significantly 


by  unsteadiness  in  highly-loaded  machines.  Note  that  (3. A. 18)  implies 


(3. A. 18a) 


r 


I ! / 
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indeed,  P(t)  = (p  “ °°)/p  = const,  when  there  is  no  unsteadiness  in  the  inlet, 
o 

(Note  also  that  if  the  only  time-dependence  is  that  associated  with 
steady  rotation  of  the  blade  surfaces,  so  that  3ct/3t  = -<o,  the  above 
conserved  quantity  becomes 


p p 

— + G (o  = — - corvQ  = const. , 
P P 0 


(3. A. 18b) 


i.e.,  it  is  then  precisely  the  "rotary  stagnation  pressure",  as  expected.) 

The  equations  determining  the  (aQ,  i|>)  surfaces  provide  a useful 
alternative  formulation  of  the  mean  flow  problem,  which  can  be  used,  for 
example,  as  a complement  to  Eqs.  (3. A. 6, 7)  for  If.  Given  that  "V  is  described 
in  the  present  case  by  (3. A. 15)  at  any  instant,  while  curl  V = = 

Va  x VG  from  (3. A. 4),  we  have 


fl  .•  Va  = 0 ; S=  «,Tl,t) 


or 


-5/r  - 5 -Vf  - i (|£  If  - If  |f>  - V6.(VGxVf) 


But  T)  = (3u/3z  - 3w/3r) , while  the  mean  meridional  velocity,  (u,  w) , 

(IM  i|i)  from  (3. A.  15).  Then 
r 3z  ’ r 3r 


A*\l>  *=  r2V0  • (VG  x (Vf  - VC)) 


(3.A.19) 


where  A*ip  = V2^  - — in  cylindrical  coordinates. 


I 
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This  is  the  appropriate  equation**’^  (for  the  incompressible  case) 
determining  the  Stokes  stream  function  In  terms  of  the  mean  vorticity 
generated  at  the  blades.  This  equation  is  equivalent  to  any  of  a variety 
of  formulations  of  "axisymmetric  throughflow  theory".  In  practical 
situations  it  is  solved  numerically  (exact  analytic  solutions  are  very  rare) ; 
for  our  present  purposes,  we  regard  its  solution  as  known  or  readily  at 
hand,  given  the  rhs  or  sufficient  information  to  determine  it:  (torque 
schedule , mean  swirl , . . . ) . 

At  the  same  time,  the  mean  swirl  surfaces,  aQ  = 9 - fQ(r,z,t),  are 
determined,  at  each  instant,  from  the  requirement  V • VaQ  = co;  that  is, 

rV  • Vf  = v„  - cor 
- o 0 


or  by  characteristic  integration  at  constant  if>(and  t) 


dz*  , -G  \ . r-o 
— ( -j  - u)  + f (t) 

w r 

const,  ip 


(3.A.20  ) 


In  order  that  the  relative  camber  displacement,  C(r,z,t)  in  (3. A. 16), 
remain  bounded  for  all  time,  even  when  the  a-surfaces  are  rotating,  we 
must  choose  f°  (t)  = cot.  Then 

(-^f  - co)  . (3.A.21) 

r 


a «*  6 - cot  + 
o 


r z 
dz' 

W 


1 
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•When  there  is  no  rotation,  of  course,  o>  = 0. 

We  now  see  explicitly,  either  from  the  formulation  in  (3. A. 19)  or  from 
that  in  (3. A. 7),  that  the  mean  flow  "reference"  problem  is  coupled  with  the 
perturbation  problem  (at  least)  through  the  displacement,  C,  relating  the 
a-  and  - surfaces.  The  function  C(r,z,t)  must  satisfy,  according  to 
(3. A. 16a) , 


DC 

Dt 


<v 


s • Va 

> O • 

bl,v 


(3.A.22) 


While  aQ  depends  only  on  mean-flow  quantities,  C(r,z,t)  can  be  determined 
only  by  solving  for  the  perturbation  velocity  v. 

As  anticipated,  these  results  show  that  any  unsteady  effects,  even 
though  they  may  be  "small"  in  some  sense,  directly  affect  the  mean-flow 
itself,  to  first  order  in  the  perturbations.  They  do  so  in  two  ways: 

1)  kinematically,  they  enter  into  the  "source"  terms  for  either  <j>(3.A.7) 

or  (3. A. 19)  as  a result  of  the  displacement  between  a-  and  a - surfaces; 

2)  dynamically,  they  have  a direct  effect  on  the  mean  "total  pressure" 
(see  (3. A. 18a))  and  hence  on  the  instantaneous  9-averaged  thermodynamic 
properties. 

In  the  incompressible  case  these  two  effects  of  time  dependence  can  be 
analyzed  separately,  as  we  have  shown  above.  For  compressible  flow,  however, 
the  dynamical  and  kinematic  features  are  intertwined  in  a very  interesting 
way.  In  either  case,  important  aspects  of  the  mean  flow  remain  undetermined 
until  the  three-dimensional  perturbation  problem  itself  is  solved.  This 
can  only  be  done  in  conjunction  with  self-consistent  perturbations  of  the 
mean-flow  problem  as  indicated  above,  and  in  the  following. 


. iii..',, 


IT 


III.  A ii.  THE  PERTURBATION  PROBLEM  FOR  UNSTEADY  FLOW 
IN  HIGHLY-LOADED  TURBOMACHINES;  UNIFORM  INLET  FLOW. 

With  uniform  inlet  flow  the  vorticity  appropriate  to  the  present  model 
is,  from  (3.3), 


while 


ft  - VH(a)  x VG 


ft  = Vet  x VG 


Hence  the  peturbation  vorticity,  = 9 - ft,  is  simply 


ft  = VS(a)  x VG 


(3.A.23) 


(3.A.2A) 


(3. A. 25) 


where  S(a)  = H(a)  - a is  the  sawtooth  function  introduced  briefly  below 
Eq.  (3.3).  S(a)  has  the  values  ±tt/B  on  the  pressure  (+)  and  suction  (-) 
sides  of  each  of  the  blades,  respectively.  Also,  from  the  point  of  view 
of  generalized  functions,  S(o)  = 0.  That  is,  operationally,  it  is  zero 
"in"  each  blade  or  wake  surface.  It  is  also  periodic,  S(ot)  = S (a  + -^^)  , 
and  has  zero  0-average  over  a blade-passage  interval. 

The  perturbation  velocity,  v,  corresponding  to  the  vorticity  in  (3. A. 25), 


v - V<J>  + SVG, 


(3.A.26) 


and,  since  div  y *■  div  y ■ 0,  we  require  div  v » 0.  Then 


As  we  have  seen  earlier,  with  V 


GVa  + S(a)VG,  and  V ■=  V<j>  - GVa,  we  have  both  for  the  complete 


Note  that  actually  these  Eqs.  are  not 


more  difficult  than  (3. A. 26).  However,  the  applicable  boundary  conditions 


-hub, 

shroud 


blade 

wake 


-hub, 

'shroud 


Swirl 

surfaces 


(perturbation) 
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vhere,  by  definition,  a = a©  + C(r,z,t),  and  is  still  to  be  determined. 

The  presumption  is,  of  course,  that  if  a perturbation  technique  is  at 

all  applicable,  then  not  only  are  the  swirl  surfaces,  a©,  determined 

entirely  from  the  mean  flow,  but  also  that  |Vc|  <<  |Va©|  ; in  fact,  we 

expect  N /|Va  | = 0 (1/Bn) , where  n > 1.  [One  should  not  be  misled  by 
° 3a 

the  presence  of  the  term  (&)  = ) in  (3. A. 28c);  e.g.,  co  = 0 if  the  blade 

forces  are  not  on  rotating  surfaces,  while,  generally,  (3. A. 21)  shows  that 


*=ct. 


9 


3a© 

so  (co  + — g-j-  ) is  actually  "small"  (more  precisely,  perturbative)  and  represents 
the  extent  to  which  a©-surfaces  are  "rippled"  by  inherent  unsteady  effects.] 

On  the  basis  of  these  equations  the  perturbation  procedure  can  be 
conceptualized  readily.  Equation  (3. A. 27b)  contains  no  singular  functions, 
but  the  rhs  still  involves  the  unknown  displacement  C(r,z,t).  Let  us  then 
define  a "zeroth-order"  mean  velocity 


and  corresponding  perturbation  mean  velocity,  v©  such  that 
v^Y-Vo-vi- $©)  - GVC  = v(4>©)  - G7C 


(3. A. 29) 


(3.A.30) 


22i- 
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( Clearly,  vQ  Is  a purely  meridional  vector,  has  no  0-dependence  and  ljs 
its  own  0-average.)  Thus,  if  the  net  perturbation,  relative  to 


V , is  "v",  then 
~o  ~ 


v = v + VQ  = V<f>o  - GVC  + v<t>  + SVG  = V^h  SVG  - GVC 


(3.  A.  31) 


while 


and 


co  = curl  v = VC  x VG  + VS  x VG 


Ti  - curl  V = Va  x VG 
~ o o o 


(a) 

(3.A.32) 

(b) 


Note  that  V , as  defined  by  (3. A. 29),  is  keyed  to  the  same  value  of  G 
o 

(cf.  (3.A.6))  as  is  the  actual  mean  flow,  V.  In  addition  to  having 
identical  ©-components  (=  -G/r) , each  is  referred  to  the  same  meridional 
vorticity : 


a -a  - (iff  > - iff)  = w x ™ 

~mo  ~m  r dz  r dr 


(3.A.33) 


However,  if  V is  to  be  useful  as  a first  approximation  (zeroth  iterate) 

for  the  mean  reference  flow,  we  must  be  able  to  use  it  to  generate  the 

a - surfaces,  to  a good  approximation,  without  prior  knowledge  of 
o 


C(r,z,t).  Thus,  we  need  also  Vq  • VaQ 


V • Va  = a>  , or 
~ o 


V • Vf  « V • Vf  - -(G)  + G/r  ) . 
mo  o ~m  o 


(3. A. 32c) 


These  requirements  are  met  adequately  if  we  choose,  as  a complement 
to  (3.A.29), 


V = - ; efl  + V0  x 
~o  r u o 


(3.A.34) 


whence  V • Vf  = - (Vf  x ) • 7\p  , and  to  fulfill  (3. A. 32c): 
~o  o o o 


fo  (r,z,t) 


z 

dz1  , __G 
_ t 2 

w r 

o 


constant  ^ 


(3.A.35) 


(compare  (3.A.20),  (3.A.21)).  Here,  wq  = - (1/r)  (3ij>o/3r)  . 

Note  that,  in  this  initial  approximation,  the  and  surfaces  do  not 
generally  coincide*.  However,  we  can  assure  that  they  do  so,  more  and  more 
as  the  perturbation  procedure  is  carried  out,  by  requiring  1 ViUb  * ^|hub 
and  *o|shroud  «■  ^js^rouj  38  boundary  conditions  in  the  equation  for  \pQ 
described  below. 

As  with  ip,  the  equation  for  is  obtained  by  requiring  consistency  with 
the  vorticity  condition,  in  this  case  (3. A. 32).  Thus,  with  • VfQ  = rio/r, 
we  find 


A*i|»  - r V6  • (VG  x Vf  ) 

o o 


(3.A.36) 


* An  exception  occurs  when  C(r,z,t)  = C(G) . 
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can  then  be  determined  from  this  equation,  subject  to  the  hub-shroud 
boundary  conditions  stated  above.  This  allows  one,  in  turn,  to  determine 
the  first  approximation  to  the  swirl  surfaces,  , through  (3. A. 35). 

We  are  now  in  a position  to  return  to  (3. A. 26a)  to  determine,  to  lowest 
order,  the  blade-to-blade  perturbations,  given  appropriate  boundary 
conditions  at  the  blades,  as  well  as  at  the  hub  and  shroud.  (The  blade 
boundary  conditions  can  be  rather  elaborate,  depending  on  the  problem  at 
hand,  but  in  any  case  must  include  (3. A. 10)  — or  its  generalization  for  a 
given  blade  thickness  profile  — as  a matter  of  consistency.)  In  treating 
the  blade-to-blade  problem,  it  is  essential,  on  the  rhs  of  (3. A. 26a)  or  at 
any  stage  (n)  of  an  iterative  or  perturbative  procedure,  to  use: 


S(ct(o)), 


(n) 

s(a  ); 


lim  a. 
n -*-  « 


(n) 


rather  than  S(a)|a=<3t  . (Otherwise,  the  perturbative  method  will  diverge.) 

n 

On  solution  of  this  problem,  C(r,z,t)  ...  or  its  corresponding  iterate, 

C^n^(r,z,t)  is  available  (see  (3. A. 28c)  with  (3. A. 26a))  and  a = otQ  + C 

(or  = a^n  ^ + C^)  can  be  determined. 

° ° 

Finally,  returning  to  (3. A. 20)  and  (3. A. 19),  improved  versions  of  the 
(ttQ,  ip)  surfaces  can  be  obtained;  thes£  are  then  to  be  used  to  find 
improved  solutions  of  (3. A. 26a),  and  so  on,  ...,.  Providing  due  care  is 
taken  to  maintain  internal  consistency,  and  the  boundary  conditions  are 
appropriately  monitored,  this  procedure  converges  extremely  rapidly.  (For 
analytic  purposes  much  of  the  above  can  be  handled  by  implicit  methods.) 


I 


There  is  a useful  physical  picture  underlying  starting  the  above 

perturbative  solution  of  the  unsteady  flow  problem  at  high  loadings  in 

the  manner  outlined  above.  By  initially  concentrating  distributed 

vorticity  of  (relatively  well-known)  strength  on  the  mean-flow  swirl 

(a  ) surfaces  in  high-deflection,  turbomachinery  flows,  one  avoids  the 
o 

difficulties,  encountered  in  earlier  studies,  of  allowing  any  shed  or 
trailing  vorticity  to  be  convected  "through"  the  regions  of  concentrated 
bound  vorticity,  i.e.  "through"  what  amount  to  solid  surfaces! 

HI.  B.  CLEBSCH-HAWTHORNE  METHOD  FOR  THIN  AIRFOIL  EXAMPLE 
In  this  section  we  offer,  very  briefly,  an  example,  taken  from  the 
classical  theory  of  incompressible  unsteady  two-dimensional  flow  past 


thin  airfoils,  for  which  the  above  analysis  makes  contact  with  well-known 
for  a simpler  but  related  problem. 

Returning  to  (3.4),  we  recall  that  a key  step  in  dissecting  that  relation 
was  the  identification  of  the  "boundary-condition"-term  (3rd  term  on  l.h.s. 
of  (3.4)  vanishes)  and,  also,  the  "body-force"  termleading  to  (3.4a). 

For  the  thin  airfoil  case  (steady  or  unsteady)  the  reference  flow  is  a 
uniform,  steady  potential  flow  Vq  = (U,0,0)  in  an  (x,y,z)  cartesian  system. 
[Conventionally,  x -*■  (streamwise)  , y (spanwise)  , z -*■  (normal  - to  airfoil)]. 
Then,  on  linearizing 


Dt  Dt  at  U 3x',_  Dt 


Again,  the  controlling  equations  are,  for  this  example. 


(3.B.1) 


_ _ v X n =:  -V  (— ) + F ; div  V = 0 
9t  p ~ 


In  this  case,  VpQ  = 0,  except  for  perturbations.  The  flow  is,  correspond 
ingly,  potential  (curl  V = 0)  aside  from  perturbations  introduced  by  the 
airfoil  and  in  its  wake. 

Thus,  if  V = V<J>  + A as  before,  and  if  the  vorticity  is  written 


n = VH(z)  X VG  = 6 (z) (e  x VG) 
~ z 


(3.B.3) 


we  find  S2  to  be  concentrated  on  the  plane  z = 0,  and  "G”  is,  itself,  a 
perturbation  quantity  — still  to  be  determined.  This  is  precisely  the 
picture  offered  in  thin  airfoil  theory. 

Making  the  same  expansions  as  in  section  III,  we  have 


DH , . 
(z) 

Dt 


VG  - ^ VH(z) 


-7(ir + It + H !?> + 


(3.B.4) 


But,  from  the  (exact)  boundary  condition,  the  first  term  vanishes. 
Identifying  F(x,t)  as  before,  viz., 


v — PG  f . _ DG  , < . 
F = - — VH(z)  = - ^ 6(z) 


ni- 


a a D G 

- + «(*>  5 - «(z) 


WT 


i 


*o  , 31  , v 3Gv  . 
(T  + 3t  + H(z)  Te>  = ct- 


(3.B.6) 


3G 


as  the  "Bernoulli  constant".  (Note  the  jump  in  the  5—  term  across 

O t 

z = 0.) 

On  the  other  hand,  the  physical  meaning  of  F is,  in  this  case,  such 
that  it  must  balance  the  instantaneous  jump  in  static  pressure  ("loading") 
at  the  airfoil  surface: 


F «=  (x.t)  3(z)  = 


VH 


Thus,  from  (3.B.5) , 


D G 


o_  ^ _ 3G  3G  _ Ap(x,t) 

Dt  3t  3x  p 


or,  if  Ap(x,t)  = Apw(x)e1/i)t  and  G(x,t)  = G<J0(x)etot 


3G, . 


-JttG  - + (x) 

to  3x  p 


and,  if  the  leading  edge  of  the  airfoil  is  at  x = 0, 


G0)(x)  " - 

• 

+i  cox/u 
e 

“X 

x'n 

► * 

0 

(3.B.7) 


(3.B.8) 


(3.B.9) 


Then 
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VV*>  ■ % - sir  - ¥ *■’  0 

o 


jcoCx-x') 


and,  identifying  yfi  = YB0UND  = (Ap(x,t)/pU)  , we  have 


■ *7  «*>  -y^c-r  d*'  rm  eto(,I‘x’)/u 


This  is  the  classical  thin-airfoil  result,  distinguishing  the  bound 
vorticity  on  the  airfoil  and  the  free,  "shed"  vorticity  in  the  wake, 
x > CHORD. 

With  (3.B.10),  one  can  determine  the  relationship  between  the 
instantaneous  blade  loading,  the  vorticity  on  and  off  the  blades,  the 
corresponding  velocity  field,  and  using  (3.B.6),  the  pressure  field. 

These  results  agree  precisely  with  those  obtained  from  "classical" 
procedures  (including  for  example,  use  of  the  Biot-Savart  Law  - which  is, 
unfortunately,  extremely  awkward  in  applications  to  internal  flows  such 
as  those  typical  to  turbomachinery) . 
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FUTURE  STUDIES 

* 

As  a result  of  the  work  accomplished  so  far,  there  are  many  lines 
open  to  future  research.  Examples  are  extension  of  present  work  to  the 
case  of  the  behavior  of  blade-wakes  (consisting  of  the  Beltrami  vorticity 
as  well  as  that  from  boundary  layers  swept  downstream  of  the  blades)  in 
flows  with  the  reference  flow  departing  substantially  from  free-vortex 
swirl,  and  of  more  importance,  extensions  through  highly-loaded  blade  rows 

: 

I 

of  finite  chord  length,  c.  Comparison  of  the  present  theory  with  available 
experimental  data  (obtained  by  work  sponsored  by  the  AFOSR  in  UTRC  and  at 
M.I.T.  Blow-Down  Compressor  Facility)  are  currently  under  way.  Such 
comparisons  will  undoubtedly  help  to  further  our  understanding  of  the  three- 
dimensional  aspects  of  turbomachinery  flow,  and  aid  us  in  assessing  the 
importance  of  the  finiteness  of  the  blade-chord,  thus  leading  to  improvements 
and  extensions  of  our  analytical  work. 

Another  area  in  which  a potential  exists  for  a major  contribution  is 
in  the  combination  of  the  analytical  approach  described  above  with  the 
newly  developing  time-marching  numerical  schemes  for  studying  turbomachinery 
flows  - steady  flow  and  unsteady  flow  in  future.  In  the  numerical  solution 
of  turbomachinery  flow,  for  practical  purposes,  it  is  almost  always  the  case 
that  uniform  boundary  conditions  are  imposed  at  approximately  one  chord 
upstream  and  downstream  of  the  blade-row.  Theoretical  studies  have  shown 
that  upstream  acoustic  radiation  occurs  in  transonic  blade-rows,  and  the 
existence  of  oscillating  disturbances  downstream  of  the  blade-row. 

Thus,  the  upstream  and  downstream  boundary  conditions  are  never  actually  uniform. 


/ 
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except  in  the  truly  "incompressible"  limit.  Thereby,  it  emerges  that 
our  analytic  work  can  be  used  to  provide  the  near  upstream  and  down- 
stream boundary  conditions,  thus  allowing  more  precise  numerical  study 
of  blade-passage  regions.  Such  a numerical-analytical  "hybrid"  scheme 
helps  to  put  numerical  solutions  of  turbomachinery  flows  on  sound 
mathematical  and  scientific  ground.  Future  studies  will  certainly  include 
analysis  of  unsteady  flows  in  turbomachines  including,  now,  of  that  in 
highly  loaded  blade  rows. 

In  contrast  to  our  earlier  works  where  lifting  line  models  and  activator 
disc  models  are  used  (see  Introduction) , the  present  and  future  emphasis 
will  be  on  the  more  realistic  models  of  highly-loaded  blade  rows  of  finite 
chord  length.  It  is  expected  that  the  later  extensions  will  cover  both 
steady  and  unsteady  problems,  and  with  particular  emphasis  on  "three- 
dimensional  design"  of  highly  loaded  rotors.  As  illustrated  in  Sections  III 
and  III. A.,  some  of  these  extensions  are  already  in  progress.  In  what 
follows,  a brief  account  is  given  of  the  three-dimensional  problems  of 
steady  and  unsteady  flows  through  highly-loaded  rotor  of  finite  chord  length, 
with  possible  application  to  the  rotor  design  problem. 

In  the  context  of  the  design  problem,  we  seek  the  determination  of  a 
blade  shape,  given  a specified  loading  (i.e.,  F) . The  basic  notion  then 
involves  the  finding  of  appropriate  blade  shapes  by  distributing  vorticity 


on  the  a -surfaces  so  as  to  create  the  required  local  F.  The  specification 
o ~ 


of  F,  aside  from  allowing  the  determination  of  the  stream  surfaces  ip  and 


the  a surfaces  from  the  mean  flow  equations,  enables  one  to  compute  the 


strength  of  the  distributed  bound  vorticity  on  the  aQ  surfaces.  This  then 


permits  the  actual  camber  lines  relative  to  the  surfaces  to  be  calculated. 


4 
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The  advantages  of  such  an  approach  for  the  design  problem  are:  (i)  large 
deflections  are  included  in  the  mean,  azimuthally-averaged  flow,  and 
(ii)  any  shed  vorticity  emanating  from  the  bound  vortex  system  will  be 
carried  along  the  surfaces,  thus  eliminating  the  problem  of  earlier 
linearized  theories  of  having  such  free  vortices  crossing  the  airfoil  surfaces. 
Furthermore,  such  an  approach  makes  possible  the  interpretation  of  the  physical 
phenomena  associated  with  the  flow.  It  can  also  be  adapted  to  the  'analysis' 
problem,  and  thus  its  use  for  predictive  and  interpretive  purposes.  As 
shown  in  Section  III,  the  general  formulations  of  the  b lade- to-b lade 
equations  (or  rotor  passage  equations)  in  terms  of  the  multiple  Glebesch- 
Hawthorne  transformation  may  turn  out  to  be  an  elegant  and  practical  way 
to  study  the  unsteady  internal  flow  in  turbomachines.  Two  possible  approaches 
to  the  solution  are  available:  (i)  an  iterative  method,  and  (ii)  an  implicit 
perturbative  scheme.  This  arises  because  the  mean  reference  flow  (time- 
average  and/or  azimuthal  average)  is  coupled  with  the  perturbation  flow 
through  the  displacement,  C (a  relative  camber),  relating  the  a-  and  a 

° 

surfaces,  together  with  the  various  constraints  on  the  flow.  Section  III.B 
shows  how  the  Clebesch-Hawthorne  approach  recovers  the  results  of  the 
classical  thin-airfoil  theory  elegantly. 

Preliminary  works  on  unsteady  flows  in  turbomachines  using  the  present 
approach  has  indicated  the  significance  of  three-dimensional  effects  with 
regard  to  blade  flutter  boundaries,  even  though,  up  til  now,  only  lightly- 
loaded  rectilinear  cascades- are being  considered.  As  mentioned,  earlier, 
future  studies  will  attempt  to  treat  the  unsteady  flow  through  highly-loaded 
blade-rows,  and  to  make  possible  combined  numerical  analytical  study  of 
unsteady  flow  phenomena  in  turbomachines. 
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FIGURE  CAPTIONS 

Radial  distribution  of  flow  angle  deflection  given  by  the 
free  vortex  solution,  the  axisymmetric  activator  disc  solution 
and  the  present  three-dimensional  Beltrami  flow  solution. 

Comparison  of  the  blade-to-blade  radial  Mach  number  as  given  by 
the  Beltrami  flow  theory  with  experimental  data  obtained  at  M.I.T. 
Blow-Down  Compressor  Facility. 

Comparison  of  the  blade-to-blade  pitchwise  Mach  number  as  given 
by  the  Beltrami  flow  theory  with  the  experimental  data  obtained  at 
M.I.T.  Blow-Down  Compressor  Facility. 

Comparison  of  the  blade-to-blade  axial  Mach  number  as  given 
by  the  Beltrami  flow  theory  with  the  experimental  data  obtained  at 
M.I.T.  Blow-Down  Compressor  Facility. 

Comparison  of  the  blade-to-blade  radial  Mach  number  as  given 
by  the  Beltrami  flow  theory  with  the  experimental  data  obtained  at  „ 
M.I.T.  Blow-Down  Compressor  Facility. 

Comparison  of  the  blade-to-blade  pitchwise  Mach  number  as  given 
by  the  Beltrami  flow  theory  with  the  experimental  data  obtained  at 
M.I.T.  Blow-Down  Compressor  “Facility. 

Comparison  of  the  blade-to-blade  axial  Mach  number  as  given 
by  the  Beltrami  flow  theory  with  the  experimental  data  obtained  at 
M.I.T.  Blow-Down  Compressor  Facility. 
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Fig.  8 

Fig.  9 

Fig.  10 

Fig.  11 

Fig.  12 

Fig.  13 
Fig.  14 
Fig.  15 
Fig.  16 

Fig.  17 

Fig.  18 


Radial  distribution  of  flow  angle  deflection  given  by  the  radial 
equilibrium  solution,  the  axisymmetric  activator  disc  solution, 
the  three-dimensional  solution  with  the  circumferential  non- 
uniformity due  to  the  presence  of  wakes  with  total  pressure  defect. 
Downstream  convection  of  stagnation  pressure  distortion  at 
r/r^,  = 0.83  as  given  by  the  linearized  theory  and  experimental  data. 
Axial  velocity  distribution  at  different  axial  stations  for  radius 
r/rT  = 0.83. 

Comparison  of  downstream  static  pressure  perturbations  at  radius 
r/r^  = 0.83:  Theory  vs.  Experiment. 

Comparison  of  downstream  flow  angle  perturbations  at  radius 
r/rT  = 0.83:  Theory  vs.  Experiment. 

Downstream  variation  of  static  pressure  perturbation  at  the  tip. 
Downstream  variation  of  radial  velocity  at  radius  r/r^,  - 0.55. 
Downstream  variation  of  absolute  flow  angle  at  the  hub. 

Downstream  variation  of  static  pressure  perturbation  at  the  tip 
for  the  case  with  five  loop  type  of  stagnation  pressure  distortion 
introduced  far  upstream. 

Pitchwise  variation  of  axial  velocity  at  various  radial  stations 
as  given  by  the  strong  Beltrami  flow  theory  and  the  M.I.T. 

Blow-Down  Compressor  Facility. 

Pitchwise  variation  of  tangential  velocity  at  various  radial 
stations  as  given  by  the  strong  Beltrami  flow  theory  and  the 
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M.I.T.  Blow-Down  Compressor  Facility. 
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